Report for Model 10

— for both jointly and separately

2024-10-08

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 300000 iteration;
  • 200000 burn-in;
  • 100 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum      1.005       1.04
## LDevsum      1.002       1.03
## dh0          1.190       1.68
## dh1          1.137       1.49
## dh2          1.039       1.09
## dl0          1.195       1.72
## dl1          1.062       1.08
## dl2          1.043       1.11
## dl3          0.997       1.00
## 
## Multivariate psrf
## 
## 1.18

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1221.44081 22086.7640
DIC3 1159.61960 22148.1564
PWAIC 44.35335 240.4214
WAIC 1187.60308 22170.2218

2. Separate Modeling of High-Level

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 300000 iteration;
  • 200000 burn-in;
  • 100 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.13       1.46
## LDevsum       1.11       1.42
## dh0           3.95       8.42
## dh1          14.95      33.82
## dh2          15.70      49.60
## dl0           1.80       3.31
## dl1           1.05       1.08
## dl2           1.57       2.96
## dl3           1.28       1.89
## 
## Multivariate psrf
## 
## 11.9

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

H2
DIC 1385.96801
DIC3 1247.70411
PWAIC 90.14037
WAIC 1315.56019

3. Separate Modeling for Low-level

3.1 Trace plots for convergence check

The current MCMC setting is:

  • 300000 iteration;
  • 200000 burn-in;
  • 100 thinning.

3.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## Devsum        1.25       1.82
## LDevsum       1.11       1.42
## dl0           1.80       3.31
## dl1           1.05       1.08
## dl2           1.57       2.96
## dl3           1.28       1.89
## 
## Multivariate psrf
## 
## 1.51

3.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

3.4 WAIC results

L4
DIC 22221.4939
DIC3 22238.8382
PWAIC 298.0137
WAIC 22271.3988